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The most unusual and debated feature of high- 
temperature superconductivity (HTSC) is the pseudo- 
gap (PG) phase, the high-temperature phase in the 
underdoped region of the phase diagram between the de- 
struction of superconductivity at and the pseudogap 
temperature T*.^'^ While the zero-temperature phase 
diagram of HTSC is relatively well understood, there 
is currently much experimental and theoretical inter- 
est in the intermediate-temperature PG phase. In this 
phase, several surprising experimental features appear: 
e.g. angle-resolved photoemission spectroscopy (ARPES) 
shows a state which is partially gapped on the experimen- 
tal Fermi snri&ceM^^^l^ 

Theoretically, the low-temperature physics of HTSC 
is well described by variational wavefunctions of the 
t-J mode L^°d^d^d3, 14^15, 16 rpj^g antiferromagnetic par- 
ent state at half filling is destroyed as doping is in- 
creased. Due to the gain of spin-exchange energy in 
the Gutzwiller-projected state, a rf-wave mean-field or- 
der is favored away from half filling. The characteris- 
tic dome for the off-diagonal order can be reproduced 
variationallyiiiii^ Low-lying Gutzwiller-projected quasi- 
particle excitations reproduce well many experimental 
features^ii^^>20^>22^>2i^ The main disadvantage of 
the variational approach is that it is a zero-temperature 
theory and cannot easily be extended to finite tempera- 
ture or to high-energy excitations.^'^ 

Many years ago, it was noticed that there is a redun- 
dant description of Gutzwiller-projected fermionic wave- 
functions exactly at half filling, parametrized by local 
SU(2) rotationsi^i^ Away from half filling, this redun- 
dancy is lifted. Later, Wen and Lee et al. proposed 
a slave-boson field theory, where the redundancy is pro- 
moted to a dynamical SU(2) gauge theory away from half 
filling,^i2ai2a The advantage of the SU(2) slave-boson ap- 
proach is that it incorporates strong correlations when 
gauge fluctuations around the mean-field saddle points 
are included. Integrating over all gauge-field configura- 
tions in this approach enforces the Gutzwiller constraint 
Hi < 2. The slave-boson mean-field theory is then not 
restricted to low temperatures. 



The SU(2) approach to the t-J model predicts that 
a state with staggered magnetic fluxes through the pla- 
quettes of the square lattice should be close in energy 
to the d-wave superconductor at low doping i^i^i^!^ 
In fact, a staggered SU(2) rotation on nearest-neighbor 
sites transforms the d-wave superconductor (SC) into the 
staggered-flux (SF) state. These two states are identical 
at half flUing. At small doping, one expects the local 
symmetry to be weakly broken, and the SU(2) rotation 
provides a route to construct a low-lying nonsupercon- 
ducting variational state for the weakly doped t-J model. 
This led to the proposal by Wen and Lee that the pure SF 
state should be realized in the vortex cores of HTSCi^ 
Indeed, it was confirmed numerically that the Gutzwiller- 
projected SF state is a very competitive ground state 
of the t-J modeli^ Further support for the SU(2) ap- 
proach came from the discovery of SF correlations in the 
Gutzwiller-projected d-wave superconductor.'^'* 

In this paper, we restrict ourselves to the so-called 
"staggered 0-mode" which interpolates between the SC 
and the SF states As the temperature is increased 
through Tc in the underdoped compounds, vortices pro- 
liferate and eventually destroy the phase coherence. In 
order to form energetically inexpensive vortices in the 
superconductor, the order parameter rotates to the SF 
state inside the cores.— However, in contrast to vortex 
cores, we do not expect a pure SF state to be realized in 
the bulk. The PG state should be viewed as a thermal 
average over different intermediate states between the SF 
and the SC state, parametrized by appropriate SU(2) 
rotations. 

In the superconducting phase at low temperature, it is 
sufficient to include Gaussian fiuctuations away from the 
superconducting state. In this framework, Honerkamp 
and Lee found that coupling to the Gaussian 0-mode 
strongly depletes the antinodal quasiparticles.'^- This 
is in contrast to zero temperature, where Gutzwiller- 
projected excitations show rather weak reduction of spec- 
tral weight in the antinodal region as shown by the 
authorsi^ At temperatures between Tc and T*, strong 
fluctuations toward the SF state are expected to affect 
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the electronic spectral functions even more. In the 
present work, we are interested in the electron spectral 
intensities in the pseudogap region, i.e. in the presence 
of large fluctuations of the order parameter between the 
SC and the SF states. 

Our model bears some similarity to the cr-model ap- 
proach for the SU(2) gauge theory of the t-J model, intro- 
duced by Lcc et al}^-^'^ In contrast to these authors, we 
do not use a self-consistent mean-field treatment, but we 
consider an effective model with input from Gutzwiller- 
projected variational wavefunctions of the t-J model. 

A complementary study was conducted by Honerkamp 
and Lee who considered SU(2) fluctuations in an inho- 
mogeneous vortex liquid^^^ These authors computed the 
density of states and helicity modulus, and found that a 
dilute liquid of SF vortices would account for the large 
Nernst signal observed in the pseudogap phasci^ In the 
present paper, we are particularly interested in the impli- 
cations of the fluctuating-staggered-flux scenario for the 
ARPES spectra. 

Finally, let us note that our model concerns the 
low-energy spectra of cuprate superconductors, |w| < 
200 meV. The interesting high-energy anomalies 
(Iwl ~ 0.4 - 1 eV) which were discovered in recent 
experimental^"^ and theoretical"*^ works are not in the 
scope of the current discussion. 

This paper is organized in the following way. In 
Sec. m we introduce the model and describe the observ- 
able (spectral function) that we want to study. In Sec. [Ill 
we give a detailed account on the spectra of the pure 
(unaveraged) states. Finally, in Sec. IIIII we present our 
results on the spectral functions averaged over the order- 
parameter space and in Sec. IIVI we discuss the experi- 
mental implications. 



I. MODEL 

The local SU(2) rotation for the t-J modelii is conve- 
niently written using spinon doublets in the usual nota- 
tion, = (c|,c|). In terms of these doublets, the SF 
state is defined by the mean- field Hamiltonian Hsf = 
'^lUij^ipj - ^iJ2^^4^l<73^^ where Uf/ = -xts - 
jA(— , (Ta are the Pauli matrices, and the sum {i, j) 
is taken over pairs of nearest-neighbor sites. In this pa- 
per, we restrict ourselves to the following SU(2) rotation: 
Uij gjUijgj, with gj = exp[i(— | cri]. Note that 
for 9 — ^, the SF Hamiltonian is rotated to a d-wave 
superconductor, Uf^^ Ufj^ = -x^s -I- A(-)*-+^-(Ti. 
IjSC represents the d-wave mean field with pairing A. 
U^^ is the nonsuperconducting mean field with stag- 
gered U(l) fluxes equal to ±4arctan(A/x) through each 
plaquette of the square lattice. The intermediate states 
for general 9 contain both d-wave pairing and staggered 
fluxes. 

We now consider the mean-fleld Hamiltonian at the 



intermediate values of 9 between and tt: 
HMFi9) = J2^Plglm^/9,m] 

~X'Y1 tplc^Stp] ~t^Yl ^J^3^» • 

As usual, the chemical potential fj, is added to enforce 
the desired average particle number. We have also added 
a phenomenological next-nearest-neighbor hopping x'- 
Note that the parameters Xi x'l ^iid ^ of Hamiltonian 
([T|) are the effective parameters describing the variational 
ground state and quasiparticle spectrum of the t-J model, 
for the physically relevant value t ZJ: The hopping x 
only weakly depends on doping (at small doping) and is 
approximately given by x — '^/S — 100 meV J^'^^i'^ At 
10% doping, A decreases slightly from Ag^^/2 — 0.25x 
in the SC state to Ae=o - 0.2x in the SF state,2^i^ 

The value of the next-nearest-neighbor hopping is 
taken to be x' — ~0-3Xi to mimick the experimental 
Fermi surface observed in cuprates. Earlier studies of 
Gutzwiller-projected wavefunctions suggest that such an 
effective next-nearest-neighbor hopping may appear in 
the underdoped region as a consequence of strong corre- 
lations, even in the absence of the term in the physical 
HamiltonianfiSii^ Note that we keep this term unrotated 
in Eq. (P). 

In our model, physical quantities at finite tempera- 
ture are given by an appropriate functional integral over 
the mean- field parameters Uij , weighted by a free energy 
which is almost flat in the directions parametrized by gj . 
As indicated earlier, we restrict our study to staggered 
SU(2) rotations parametrized by the angle 9. At the 
same time, we neglect the amplitude fluctuations of A, 
since the energy scale associated with these fluctuations 
is high: of the order of T* , in our approach. On the other 
hand, the energy scale Sc = Esf — Esc, responsible for 
the 0- fluctuations, is much lower: at 10% doping it is 
estimated as Ec — 0.02J ~ 30K (per lattice site) from 
variational Monte Carlo calculationsi^ 

The free energy describing classical fluctuations of 9 
contains a 0-dependent "condensation energy" (of the or- 
der Ec) and a gradient terrri^^ p{'V9)^. We assume a sit- 
uation where the resulting correlation length ^ = \/ pj^c 
is much larger than one lattice spacing,''^'''^ In this case, 
the characteristic temperature, below which the conden- 
sation energy selects the superconducting state over the 
staggered-flux one, is ~ p/ln^ (the same scale deter- 
mines the temperature of the Kosterlitz-Thouless-type 
transition).— For temperatures above that scale, but be- 
low p, 

p/lne<T<p, (2) 

the order parameter slowly varies in space and takes all 
possible values related by the SU(2) rotation. Therefore, 
in this temperature range, we can approximate the clas- 
sical fluctuations by an equal-weight statistical average 
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over the uniform states with all possible values of 9. The 
corresponding integration measure for is d{cos9), 
inherited from the invariant measure on SU(2). 

We calculate the spectral function ^^(w) — 
— ilmGfc(a; + iT) where G is the single-particle Green's 
function^ of Hmf, Eq. (P). Note that Lu is measured 
with respect to the Fermi energy throughout this paper. 
As explained above, the spectra in the pseudogap phase 
are modeled by the averages of this spectral function over 
the order-parameter space, Ak{uj) = d{coa9)A^{uj). 
The spectral functions of the pure states are sums 

of delta functions. After averaging over 9, the spectral 
functions acquire an intrinsic width. In addition, we in- 
troduce a lifetime broadening T to make the figures more 
readable. 



state. However, the points on the superconducting Fermi 
surface where it crosses the diagonal (tt, 0)-(0, tt) do not 
move as 9 is changed. We will call them SU(2) points^ 
because at these points, the full SU(2) symmetry is in- 
tact even away from half filling. We will comment more 
on this later. As we decrease 0, the SC gap, symmetric 
with respect to the Fermi level, decreases and closes at 
9 = [Asc = 2(cosfca; — cosfcy) A sin0]. At the same 
time, an SF gap opens on the diagonal (tt, 0)-(0, tt) at 
the energy uj ~ —jl [we define fi = ^ — 2%' cos kx cos ky\. 
The SF gap value is I^sf — 2{coskx — cos ky) A cos 9. 
The spectral weight is transferred among the four bands 
and all of them gain intensity in the intermediate states. 
However, in most parts of the zone, there is only a single 
strong band. 




4 




3jr 
4 



3jL 
4 





(c), ■ 































In 

4 



\ 



\ 



\ 



A 



I 







3>r 
4 







4 



FIG. 1: Contour plot of the spectral function at the Fermi 
energy of the pure states, A^{0). Doping is 10% and we use 
a lifetime broadening F = 0.2x- The dashed line represents 
the location of the Fermi surface (where the Green's function 
changes sign^). From (a) to (d) we have cos 6 = 0, 1/3, 2/3, 1. 
Plot (a) shows the d-wave superconductor, (d) displays the 
pure staggered-flux state. 



II. PURE STATES 

In order to understand the averaged spectral function, 
we first outline how the intermediate states evolve as 9 is 
decreased from ^ = f (SC state) to 6* = (SF state). 
In Fig. [H we plot the Fermi surface (more precisely 
the Luttinger surface^) and the spectral function at the 
Fermi energy, A^{0). As the parameter 9 is decreased 
from ^ , the superconducting Fermi surface gradually de- 
forms to the well-known pocket around (f , f ) of the SF 




FIG. 2: Scheme of the different regions of the Fermi surface. 
The SC and SF gaps open near the node (N point; (p = 0) 
where they are small and do not overlap. The Fermi surface 
appears as a gap less arc in region I. In region II, the two gaps 
start to overlap and form an effective gap which is shifted 
upward in energy (vertical arrows). The effective gap comes 
down in energy as we move toward the antinode in region II. 
Exactly at the SU(2) points on the diagonal (0, 7r)-(7r, 0), the 
effective gap is symmetric. Beyond the SU(2) points (re- 
gion III), the midgap is shifted below the Fermi energy. 
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The SU(2) points mentioned before belong in fact to 
SU(2) surfaces (see illustration in Fig. [J) where jl ~ 
ji — 2x' cos kx cos ky = 0. On these surfaces, the full 
SU(2) symmetry is intact even away from half filling, in 
the sense that the mean-field spectra are degenerate and 
independent of (if we neglect the weak dependency of 
A on 0; see Sec. |T] and Remark ^43). 

A schematic plot of the band structure and an illustra- 
tion of the spectral- weight transfer as we go from the SC 
to the SF state is shown in Figs.[3l|4l and [5] on cuts par- 
allel to the diagonal (0, 0)-(7r, tt). The behavior is quali- 
tatively similar for all parallel cuts. The strong weights 
stay on the respective bands as they continuously move, 
except in a small stripe between the diagonal (0, 7r)-(7r, 0) 
and the superconducting Fermi surface, outside the SF 
pocket (regions II and III in Fig. [5]). In region II (// < 0), 
the strong SC band at positive energy transfers some of 
its weight to the SF band at negative energy (see Fig.[4|). 
Here, the midpoint of the SF bands lies at positive en- 
ergy. In region III (/i > 0), the strong SC band at nega- 
tive energy transfers its weight to the SF band at positive 
energy (see Fig.[S]). The midpoint of the SF bands is now 
shifted below the Fermi energy. 



* * * 



FIG. 3: Schematic evolution of the spectra along a cut parallel 
to the nodal direction, inside the pocket (through region I in 
Fig. [21 e.g. cut b in Fig. |6]). The dot size is proportional to 
the spectral weight. For upper left, upper right, lower left, 
lower right we have cosS = 0,1/3,2/3,1. Upper left is the 
superconducting state, lower right is the staggered-flux state. 
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FIG. 5: Same plot as in Fig. O but for a cut outside the 
SU(2) point (through region III in Fig.[2l e.g. cut d in Fig. [6}. 



III. AVERAGED STATE 

The gap considerations in Sec. [Til help now to under- 
stand the spectral properties of the averaged pseudogap 
state. If we neglect the weak dependency A(^), the 
average gaps in the pocket region (region I in Fig. [2]) 
may be estimated as {Asc) — f A(cosfca; — cos ky) and 
(Asf) — A(cosfca; — cosky). In the region outside the 
pocket, the midgap energy of an effective gap is approx- 
imately given by —fl/2. On the other hand, if we are 
strict with the definition of the effective gap and only 
consider the truly excitation-free region, then we come 
to a picture with a gapless arc in region I and an opening 
of an effective gap when the two gaps start to overlap in 
region II. This effective gap opens above the Fermi en- 
ergy and comes down as we move toward (tt, 0). At the 
SU(2) point, it is symmetric around the Fermi energy. 
Moving further out in region III, we find an effective gap 
with midgap below the Fermi energy (see Fig.[2|). 

In Fig. [5] we plot the averaged spectral function at the 
Fermi energy, ^^(0), with a quasiparticle lifetime broad- 
ening of F = 0.2%. This spectral intensity resembles the 
superconducting one, in contrast to the pure SF pocket 
[see Fig.[Tfd)]. In fact, the absence of pocket-like features 
( "turn in" of the arc) in ARPES measurements was used 
in Ref. 54 as an argument against the (static) staggered- 
flux state. However, in our averaged state this "turn in" 
is completely washed out by the fluctuations toward the 
SC state, and the Fermi arc closely follows the SC Fermi 
surface, consistent with ARPES data. 
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FIG. 4: Same plot as in Fig. [3] but for a cut outside the 
pocket (through region II in Fig. [5J e.g. cut c in Fig. [S]). 
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FIG. 6: Averaged spectral function at the Fermi energy, 
Ak{0). Doping is 10% and we use a lifetime broadening 
r = 0.2x- The dashed lines are the Fermi surfaces of the 
SF and SC states, respectively. The full spectra on the cuts 
a to d are given in Figs. I7I10I 




FIG. 7: (Color online) Averaged spectral function along a 
cut parallel to the nodal direction in the Brillouin zone, cut 
a in Fig. [6l The spectra are set off in y direction by kx 
(in units of tt). kx goes from 0.25 (lowest curve) to 0.55 
(highest curve) in steps of 0.01. The curve at approxi- 
mately hp is indicated on the right (red online). We use 
a lifetime broadening F = 0.12x. The parameters used are 
A{e) = ^(0.2 cos 61)2 + (0.25 sin 6')2 and doping is 10%. The 
energy to is given in units of 2x — 200 meV, the intensities are 
in arbitrary units. 



FIG. 8: (Color online) Same plot as in Fig. [7] but on cut b in 
Fig-El kx goes from 0.15 to 0.45. 




FIG. 9: (Color online) Same plot as in Fig. [7] but on cut c in 
Fig. |B] kx goes from 0.05 to 0.35. An asymmetric effective 
gap with midgap above the Fermi energy is formed. 




FIG. 10: (Color online) Same plot as in Fig. [7] but on cut d 
in Fig. [6] kx goes from to 0.3. An asymmetric effective gap 
with midgap below the Fermi energy is formed. 
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The averaged spectra on the cuts a to d in Fig. [6] are 
shown in Figs. [TlfTPl In addition to the intrinsic width of 
the averaged state, we have chosen a lifetime broadening 
of r = 0.12x in these plots. From the averaged spectral 
function we can confirm what was already anticipated 
from the pure states: 

(i) In the region near the node, one can see a small 
symmetric suppression of intensity coming from 
the superconducting gap centered at the Fermi en- 
ergy and a very small suppression coming from the 
staggered- flux gap centered above the Fermi energy. 
These "gaps" are easily washed out by broadening 
effects (see Figs. [7] and H]) and may give rise to a 
Fermi arc. 

(ii) Outside the arc, the (pseudo-)gap opens asymmet- 
rically, with midgap first above the Fermi energy 
(Fig. [5]). Closer to the Brillouin-zone boundary, 
as we cross the SU(2) point, the gap becomes 
asymmetric with midgap below the Fermi energy 
(Fig. [TU]). At the SU(2) point, the gap is exactly 
symmetric (see illustration in Fig. [5]). 

(iii) The backbending spectra at the edges of the two 
gaps lead to a doubling of the bands in some loca- 
tions of the Brillouin zone (see Figs. [5] and [S]). This 
band doubling only happens for weak branches and 
at positive energy. 

Finally, let us emphasize that the asymmetry we find 
here is in the location of the two pseudogap coherence 
peaks with respect to the Fermi energy. A different asym- 
metry in the renormalization of the coherent spectral 
weights in the superconducting state at low doping has 
been reported in recent variational Monte Carlo calcula- 
tions, where the Gutzwiller constraint < 2 is taken into 
account exactlyJ^s^i^ We expect that such a spectral- 
weight asymmetry is also present in our model (if one 



includes the Gutzwiller projection), but a confirmation 
would require an extensive numerical work. 



IV. EXPERIMENTAL IMPLICATIONS 

The most striking prediction of our model, the for- 
mation of a staggered-flux gap above the Fermi energy, 
is difficult to verify directly in ARPES experiments, be- 
cause this effect only appears at positive energy, around 
uj ~ 100 meV. On the other hand, our more sub- 
tle prediction, the combination of superconducting and 
staggered-flux gaps into a single asymmetric gap, appear- 
ing in the antinodal region of the cuprates may well be 
within current experimental reach. However, it is clear 
that the widely applied energy-symmetrization of the 
photoemission intensities-'- inevitably destroys all such 
signs in the ARPES spectra. A careful explicit removal 
of temperature- and device-dependent factors from the 
ARPES intensity will be extremely important in order 
to detect these effects. We hope that our work will stim- 
ulate experimental and theoretical effort in this direction. 

Note added. After submission of this work, a possibly 
relevant ARPES study of the pseudogap phase of under- 
doped Bi2212 was published. These authors found a 
gap asymmetry in the experimental spectral function 
which is in agreement with our findings. An alterna- 
tive approach to ours was able to fit these data with a 
phenomenological spectral functioui^ 
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